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We consider the directed percolation process as a prototype of a system displaying a non- 
equilibrium phase transition into an absorbing state. The model is in a critical state when the 
activation probability is adjusted at some precise value p c . It is known that this criticality is lost as 
soon as the probability to activate sites for the first time p\ is changed. We show here that critical- 
ity can be restored by "compensating" the change in pi by a change of the second time activation 
probability pi in the opposite direction. At compensation, we observe that the bulk exponent of 
the process coincide with those of a normal directed percolation process. Instead, the spreading 
exponents are changed, and take values that depend continuously on the pair p\-p2- We interpret 
this situation by recognizing that the model with modified initial probabilities is one with an infinite 
number of absorbing states. We mention the possibility to apply this modified directed percolation 
model to the understanding of the size distribution of seismic events. 

PACS numbers: 



I. INTRODUCTION 

There are many systems in nature that upon a contin- 
uous input of energy, react by sudden releases of the ac- 
cumulated energy in the form of discrete events, that we 
call avalanches in general. Examples are the dynamics of 
sand piles, magnetic domain inversions in ferromagnets, 
stress release on the earth crust in the form of earth- 
quakes, and many others. A remarkable characteristic of 
most of these realizations is the fact that the size distri- 
bution of the avalanches may display power laws that are 
a manifestation of the lack on an intrinsic spatial scale 
in the system. 

The theoretical analysis of such a variety of different 
processes has focused in the common features of the prob- 
lem, and has tried to isolate the minimum necessary in- 
gredients to explain the phenomenology that is common 
to most realizations. There are a number of models that 
display critical behavior, and thus a power law avalanche 
size distribution. In most cases the obtained values of 
the exponents characterizing the avalanches are limited 
to a few possible values, corresponding to different uni- 
versality classes. 

One of the reference models that are studied in this 
context is Directed Percolation (DP). DP is the paradig- 
matic example of dynamical phase transitions into ab- 
sorbing states. See [TH1] for reviews. It provides an ex- 
ample of a very robust universality class with well studied 
critical behavior where power law distributed avalanches 
are generated. One of its most remarkable characteristics 
is its robustness, that makes that many different partic- 
ular models can effectively be described within the DP 
scenario. 

It has been shown that the critical properties of the 
DP transition are lost if the probability to activate a site 



for the first time is reduced with respect to the following 
probabilities^!)]. In this paper we show that in this case, 
criticality can be restored by an appropriate increase of 
some of the following probabilities, in a process that we 
call "compensation" . Several critical exponents found at 
compensation do not coincide with those of pure DP. In 
particular, a time-reversal symmetry known to be valid 
for DP is violated at compensation. Other exponents 
conserve the values they have for DP. The values of the 
exponents that are seen to change depend in a continu- 
ous way on the location of the compensation point in a 
parameter space. 

The rest of the paper is organized as follows. In Sec. 
[EH we review the critical properties of Directed Percola- 
tion, introducing the critical exponents and the scaling 
relations. 



In Sec. Ill 



we recall the results known for a 
modified first infection model and present a variant: the 
modified first attempt model. In Sec. |IV| we present our 
results about the possibility of compensation. These re- 
sults are discussed and compared with related models in 
Sec. [Vl Finally our conclusions are in Sec. |VI| 



II. 



CRITICAL BEHAVIOR OF DIRECTED 
PERCOLATION 



DP is a dynamical model defined on a lattice, where 
each site is associated with a state (active or inactive) 
which evolves in time. Two commonly considered vari- 
ants of this model are: site DP and bond DP. In site 
DP, a site on the lattice will be active at time t + 1 with 
probability p if at least one of its neighbors is active at 
time t. In bond DP, a site will be active at time t+1 
with probability 1 — (1 — p) , k being the number of its 
active neighbors at time t. The configuration with no 
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active sites is called an absorbing state because once it is 
reached, the dynamics stops. In DP the absorbing state 
is unique. 

For p small, the system is trapped in the absorbing 
state exponentially fast, while for large p, the system has 
a finite probability to remain active indefinitely. There 
exists a threshold p c at which the system is critical, 
that is its fluctuations have no upper limit. Around the 
threshold p c the system displays a non equilibrium phase 
transition from a fluctuating phase to an absorbing state. 
As for standard equilibrium phase transitions, univer- 
sal behavior and critical exponents are expected. It was 
found that both site and bond DP belong to the same 
universality class. Here we focus on bond DP on a two 
dimensional square lattice, for which p c ~ 0.287338 [7]. 

Being p the control parameter of the transition, we de- 
note the distance from criticality as A = \p — p c \. Two 
different order parameters can be defined: When the ini- 
tial condition corresponds to a fully active lattice the 
relevant question is to determine the density of active 
sites when t — > oo (the stationary state) , namely p s t . For 
p < p c , p st = 0, for p > p c , p st = A^. When at time 
t = a single site, located at r = 0, is active, a cluster 
of active sites spreads from the origin. Here the relevant 
question is to determine the probability to remain out 
of the absorbing state when t — > oo, namely Q st - For 
p < Pc, Qst = 0, for p > pc, Q st = A 13 ' . 

Similarly to the case of equilibrium phase transitions, 
when approaching criticality, a diverging length describes 
the spatial correlations £j_ ~ A~ v± , moreover in dynami- 
cal phase transition there is a characteristic scale for time 
correlations, £n ~ A _I/ n . These scales are independent of 
the observable and thus of the initial condition, while one 
expects the two distinct order parameters p s t and Q s t to 
be characterized by different exponents (3 and j3' [S]. We 
will see that other quantities display power law behavior 
with different critical exponents, however it is possible to 
write scaling relations that constrain the set of critical 
exponents to only four independent quantities. 

In practice, in numerical simulations it is convenient 
to start from the single seed initial condition and let the 
cluster evolve up to time t. To characterize the growth of 
spreading clusters, one measures the survival probability 
Q(t) and the average number of active sites at time t, 
N(t). These two quantities obey the scaling forms: 



Q(t) cx r s 9l (t/^ 
N(t) cx f'g 2 (t/a) 



(1) 

(2) 



where g\ and 172 are 1 at t = 0, gi(t > £y) — below 
threshold. When we consider surviving clusters only, we 
can measure average spatial extension of the cluster at 
time t, L d {t), and the average density p(t) of active sites 
at time t inside this region. These two quantities obey 
the scaling forms: 



p(t) cx r e g 3 {t/^) 
L(t) cx tig 4 {t/M 



(3) 
(4) 



where 53 and (74 behave similarly to g± and g 2 below 
threshold. 

Above threshold, both Q(t) and p(t) reach their respec- 
tive stationary state, Q s t and p s t, at the characteristic 
time £11, so that two scaling relations can be written: 



(5) 
(6) 



At the critical point the scale invariance predicts that if 
time is rescaled by a factor b, space should be rescaled 
by a factor b l,± ' v ^. Thus the size of a cluster grows as 
L(t) ~ t v± '"f> and a third scaling relations can be written: 



(7) 



Finally a generalized hyperscaling relation [TD] valid below 
the upper critical dimension [5] relates the four quantities 
previously defined. Namely N(t) can be expressed as the 
the sum of two contributions: the active sites of surviving 
clusters (~ p(t)L d (t)) which have probability Q(t), and 
the contribution of dead clusters. This writes as: 



N(t) = L d (t)p(t)-Q(t)+0-{l 
n = --6-5. 



Q{t)) 



(8) 



Below threshold, each cluster can be identified with an 
avalanche and dies in a finite time T. We define the size 
S of an avalanche as the total number of activations that 
occurred, and are mainly interested in its statistics, P(S), 
which is expected to follow a power law at criticality: 
P(S) ~ S~ T . The characteristic size of an avalanche is 
related to T through 



S(T) 



Q(t) 



dt~T 



(9) 



Assuming that fluctuations around this characteristic 
value are small, we can write P(S) dS ~ —Q'(T) dt where 
-Q'(T) - T- 5 - 1 stands for the rate of death. Com- 
bining the latter relation with Eq.((9]) we have P(S) ~ 

1-1- 4-A f l + y + 25 \ 

ju+n+o s \ i+>7+« J 5 a scaling relation for the expo- 
nent r can thus be written: 



1 + 77 + 25 



1 



rj- 



(10) 



The exponents and relations that wc introduced here 
are general features of all absorbing phase transitions, 
which are characterized by only four independent expo- 
nents: d, 9, z and uu . 

The exponents ft, 9, z/y , v±_ and z are called "bulk ex- 
ponents" because they can be measured both from the 
fully active initial condition, and from the single seed 
initial condition with averages performed over surviving 
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space 
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FIG. 1: 1-dimensional bond DP. Normal direction of time is 
downwards. The arrows are given once and for all and are the 
same for both panels. Final time is t = 12. Left: DP starting 
with a fully active lattice: M = 7 occupied seeds (green), 
with final density pit) = I. Right: DP with time reversed. 
In light blue, the paths which die before the end. In red, the 
paths that survive until t. There are exactly m = 3 seeds that 
participate in surviving walks: Q(t) = I 
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runs exclusively. The exponents j3',8,r] and r are called 
"spreading exponents" because they can be measured 
starting from a single seed, with averages performed over 
all runs. 

DP has an additional symmetry associated with time 
reversal, which implies that 9 = 6 [H QT]. This is 
schematically proved in Fig[T] where arrows, drawn with 
probability p, connect neighboring sites. An activated 
site at the start of an arrow activates the site at the end 
of the arrow. The key observation is that the direction 
of time is arbitrary: starting from the top is equivalent 
to starting from the bottom with reversed arrows. This 
exact relation thus reads: 

Q(t) = p(t) (11) 
5=6. (12) 

A necessary condition for this time-reversal symmetry 
is the uniqueness of the absorbing state. In a process 
with multiple absorbing states or aging, one cannot freely 
reverse the arrow of time. 

As a reference we recall the DP exponents precisely 
measured in numerical simulations [7]: 

5 = 6 = 0.4505 ± 0.005 z = 1.7660 ± 0.0016 
v\\ = 1.1325 ± 0.006 77 = 0.2295 ± 0.005. (13) 

III. FIRST INFECTION AND FIRST ATTEMPT 
MODELS 

A generalization of the DP process is the modified first 
Infection Model (IM) [3 [HI [T3] . In this variant, the 
probability to activate a site for the first time is given 
a value p\ different than the value in the following acti- 
vations (that we call p). This has been considered as a 
model to describe epidemic processes with (partial) im- 



FIG. 2: Phase diagram of the model with a first infection 
probability pi different from the subsequent reinfection prob- 
ability p su bs ■ 

munization. In this context, the activation of a site is 
called infection, and it is understood that the possibility 
of the sw&sequent reinfection probability p su bs can differ 
from the first infection probability pi due to "immuniza- 
tion" effects. The phase diagram of this problem in d = 2 
was established in [5] and is reproduced in Fig. [2j for ref- 
erence. 

DP critical behavior occurs at pi = p su h s — p c - At 
this point, Q(t), N(t), L 2 (t) and p(t) have the power 
law distribution corresponding to pure DP. In the line 
terminated in the blue points, the system experiences a 
phase transition corresponding to the so called General 
Epidemic Process (GEP). The fixed point of (bond) GEP 
is located at pi = 1/2, p = and corresponds exactly to 
the problem of bond Isotropic Percolation. 

Along the AB line, except for the unstable DP fixed 
point, the system is not critical. In particular, the sur- 
viving probability Q(t) and the size distribution of the 
avalanches P(S) acquire a stretched exponential charac- 
ter. The instability of the DP fixed point was shown in 
[5] : the renormalization flow takes one from any point in 
AB to either A or B. 

Instead of the case in which there are different proba- 
bilities for the first infections, we will focus on the case 
in which different probabilities occur for successive at- 
tempts, namely irrespective if the activation of the site 
actually occurred or not. The state is defined by the 
number of trials of activation, not the number of infec- 
tions. The reason to study this variant is that it may 
be useful to understand the avalanche size distribution 
in some models of seismic phenomena (the connection of 
DP with this problem will be discussed elsewhere). We 
will refer to the two possibilities as the Infection Model 
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(IM), and the Attempt Model (AM). The AM is a sort 
of milder modification of the original DP problem, com- 
pared to the IM. We expect the phase diagram of the AM 
to be qualitatively similar to that of the IM. 

In particular, the DP fixed point is clearly located at 
the same position, while the GEP point is slightly dif- 
ferent. For the IM it corresponds to 2-dimensional bond 
Isotropic Percolation (px = 0.5, p = 0), and for the AM 
it corresponds to 2-dimensional site Isotropic Percolation 
(pi ~ 0.592746, p — 0). Indeed, we observe that for the 
AM, when p = 0, a site can be activated only at the 
very first attempt, with probability p\ (no matter if we 
consider site or bond DP), thus the sites that are acti- 
vated once with this rule are exactly the sites activated 
in cZ-dimensional site Isotropic Percolation. 

The main difference between AM and IM is that the 
AM has a non-singular limiting behavior as p\ — > 0, lead- 
ing in particular to a finite mean event size (S) in this 
limit, whereas for the IM (S) goes to as p\ — > 0. 



IV. RECOVERING CRITIC ALITY WITH 
COMPENSATION: MODEL AND RESULTS 

The main point addressed in the present paper is to 
show that the lack of criticality generated by a value of 
the first attempt p\ smaller than p c (resp. larger) can 
be "compensated" by a larger (resp. smaller) second at- 
tempt probability p2- We will present strong numerical 
evidence that this compensation can occur, restoring crit- 
ical behavior in the system. 

In addition, a remarkable result is that at compensa- 
tion, several critical exponents of the problem, in partic- 
ular the bulk exponents 8, z and vn , take their normal DP 
values, while the spreading exponents (5, 77, r) depend on 
the precise values of pi and pi . 

We will not discuss the possibility of compensation in 
the IM since we cannot be conclusive at present. Al- 
though it seems that compensation can be obtained, nu- 
merical evidence is not enough for a discussion on the 
variation or not of the obtained critical exponents. We 
prefer to concentrate on the Attempt Model, where we 
are much more confident with the numerical results. 

We consider the case in which the first two attempts p\ 
and P2 differ from the sw&sequent ones, that from now on 
we consider to be equal to the critical DP value: Pi>2 = 
Psubs = Pc = 0.287338. 

A heuristic argument suggesting that such a compen- 
sation can result in criticality is the following. As a per- 
turbation, the relevant character of a change in p\ was 
demonstrated in [5] for the IM. The analysis presented 
there indicates that a change in pi generates qualitatively 
the same kind of perturbation (to leading order) than a 
different p\. Therefore, it is not surprising that there are 
particular combinations of p\ and P2 at which the lead- 
ing term of both perturbations cancel each other. These 
particular combinations will be the compensating pairs of 
values (pi,£>2)- However, the fact that we do not recover 
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FIG. 3: Q(t) for different choices of (pi,P2). Circles rep- 
resent the pure DP at pi = p2 — p BU b B ~ Pc- Averages are 
performed over 1.2 x 10 6 samples. Triangles represent the 
AM with P2 = 0,p su bs = Pc- From top to bottom, we used 
pi = 0.494,0.489,0.485. For pi = 0.489, Q(t) displays a clear 
power-law with S = 0.212 ± 0.007. Averages were performed 
over about 10° samples. 



the pure DP exponents indicates that higher order terms 
do not vanish, but result in a marginal perturbation. 

We present first the numerical evidence of the compen- 
sation effect. In all simulations, we started from a single 
active site (seed) a time t = that was in a state of be- 
ing attempted twice, and let the clusters grow until their 
natural death, or time t = 10 5 . To be very precise about 
our choices: a site that has been successfully infected at 
the first attempt is still in a state of being attempted just 
once. 

We investigated two pairs of compensating points 
(Pi)P2) and compare with usual DP (j>\ = P2 — Psubs = 
p c ). For the first one, we set P2 — and varied pi in or- 
der to find the critical point. In Fig. [3j we show that we 
recover the critical character of the surviving probability 
at the point ( Pl = 0.4890 ± 0.0005, p 2 = 0). 

The critical exponent 5 measured in the two cases is 
different. For pure DP we recover the DP exponent while 
for the AM point we have S = 0.212 ± 0.007. For the 
second point we set p\ = 0.01 and varied p2- In Fig. [4j 
we show the critical character of the point p\ = 0.01, P2 = 
0.6000 ± 0.0005. Similarly we find a new value for 5: 
0.535 ± 0.005. 

Let us present the critical behavior of the quantities 
related to the bulk exponents, 9,z,v\\. L(t) corresponds 
to the mean cluster width averaged over runs that survive 
until time t. In Fig. [5] we compare our data at two 
compensation points and at the DP point: we notice that 
the z exponent does not change, unlike the coefficient 
before the power law. In Fig. [6j p(t) corresponds to the 
mean density averaged over runs that survived until t. 
The density of a single run is measured as the ratio of 
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FIG. 4: Q(t) for different choices of (pi,p2). Circles represent 
the pure DP at p\ — p2 = p su bs = Pc- Squares represent the 
AM with pi — 0.01,p su bs = Pc- From top to bottom, we used 
p 2 = 0.62,0.60,0.58. For Pl = 0.600, Q(t) displays a clear 
power-law with 8 — 0.535 ± 0.005. Averages are performed 
over about 10 s samples. 



FIG. 6: p(t) for the AM at criticality. Circles represent the 
pure DP. Squares represent the AM with pi = 0.01, P2 = 
0.600, p S ubs = Pc- Triangles represent the AM with p\ = 
0.489, P2 = 0,p su bs = Pc- Dashed line corresponds to the 
exponent measured in pure DP (131. Averages are performed 
over 10 s — 10 6 samples. 
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FIG. 5: L(t) for the AM at criticality. Circles represent 
the pure DP at pi — P2 ~ Psubs = Pc- Squares represent 
the AM with pi — 0.01, p2 = 0.600, p S ubs = Pc- Triangles 
represent the AM withpi = 0.489, p2 = 0,p su bs = Pc- Dashed 
line corresponds to the exponent measured in pure DP ( ]13[ l. 
Averages are performed over 10 5 - 10 6 samples. 



the number of active sites at t over the number of sites 
that were activated at least once. Again, one may notice 
in Fig. [6]that the exponent remains unchanged between 
the different critical points. 

We want to check if v\\ changes with p\ and p 2 - To 
do this we set (pi — 0.01, P2 — 0.600) and use different 
values of p S ubs < Pc , thus varying A, and observe the 
deviation from power-law behavior in Fig. [7]ja). We con- 
sider the scaling law Eq. 0, using the value of S = 0.535 



extracted from Fig. [4] and the DP value given in Eq.( 13 1, 
we obtain a perfect collapse for the survival probability. 
This shows that i/u does not change between compensa- 
tion and DP. 



The scenario is different for the spreading exponents 
S, r\ and r. We already saw that S changes at compensa- 
tion. In addition, in Fig. [HJ the number of active sites 
averaged over all runs, N(t), is seen to depend on the 
compensation pair (pi,P2)- For the compensated point 
(pi = 0.489, pi — 0) we measure rj = 0.46 ±0.005 and for 
{pi = 0.01, p 2 = 0.600) we measure r/ = 0.152±0.005. At 
compensation, we expect the hyperscaling relation ^ to 
hold. As z and 9 are found to be constant, the only way 
to preserve this relation is to have <5+?y = d/z — 9 = const. 
This constant is 0.6800 ± 0.00002, if we refer to [7]. 
For the point [p\ — 0.489, pi — 0,p su b s = Pc)-, we find 
that S + 7] = 0.672 ± 0.012. For the other compensa- 
tion point {pi = 0.01, P2 = 0.600, p su bs = Pc), we find 
S + rj = 0.687 ± 0.01. These results are consistent with 
the expected value, for both compensation points. 

In Fig. |9j we present the probability density function 
P(S). The scaling relation (10) holds for the compen- 
sation process. In particular for the first compensation 
point (pi = 0.489, p 2 = 0), from S = 0.212 ± 0.007, we 
expect r = 1.126±0.0042 and measure r = 1.129±0.005. 
For the other compensation point we expect t = 1.318 ± 
0.003 and measure r = 1.320±0.005. These results are all 
consistent with the expected values, within our numer- 
ical precision. We see that the relations derived in the 
first section are still valid, except for the time-reversal 
symmetry which is violated, since 5 =/= 8. 
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FIG. 7: Fig. (a): Q(t) for the AM for pi = 0.01, p 2 = 0.600 
and different p su b s - From top to bottom, p su b s = 0.287338, 
0.28733, 0.28732, 0.2873, 0.2872, 0.287. Averages are per- 
formed over 8 10 6 samples. Fig. (b): We collapse this data, 
plotting Q(t) ■ t s against t/A _1/ u We used the 6 = 0.535 mea- 
sured in figure [4] and the DP value given in 1 13 1. 



FIG. 8: N(t) for the AM at criticality. Triangles represent 
the AM with p\ — 0.489, P2 = 0,p su b s = Pc- There we mea- 
sure r\ = 0.46 ± 0.005 (dashed line). Circles represent the 
pure DP at p\ — p2 — p su bs = Pc- Dashed line corresponds 
to the exponent measured in pure DP ( 13 1 Squares represent 
the AM with p\ = 0.01, P2 = 0.600, p BU bs = Pc- We measure 
77 = 0.152 ±0.005 (dashed line). Averages are performed over 
10 5 - 10 6 samples. 
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V. DISCUSSION AND RELATED MODELS 
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A generic description of the behavior of the model can 
be presented in the (pi,P2) parameter plane (Fig. 10 1. 
In this plane there is a line along which the behavior of 
the system is critical. This line passes through the DP 
point (pi = p c ,P2 = Pc)- The values of the bulk ex- 
ponents z, 9 and i/y are constant all along the line. The 
three spreading exponents <5, r\ and r change continuously 
when we move along the line, but always respect the rela- 
tions ^ and (10 1, so that there is only one independent 
exponent that changes. 5 passes from lower-than-DP val- 
ues when pi > p c > P2, to larger-than-DP values when 
Pi < Pc < P2- Out of this line, there is in general a 
stretched exponential contribution to the distribution of 
the relevant quantities of the problem. 



FIG. 9: P(S) for the AM at criticality. Circles represent 
the pure DP at pi — P2 = Psubs = Pc- We check that r = 
1.268±0.005. Squares represent the AM with pi = 0.01, p 2 = 
0.600, fsubs = Pc- We measure r = 1.32 ± 0.005. Triangles 
represent the AM with pi = 0.489, P2 = 0, p su b s = Pc- There 
we measure r = 1.126 ± 0.005. Averages are performed over 
10 5 - 10 6 samples. 



We can understand the behavior of the bulk exponents 
if we think that these exponents can be measured starting 
from a fully active lattice. In this case the evolution of 
the system coincides with that of pure DP after exactly 
two time steps. However, bulk exponents can also be 
measured on the surviving runs started from a single ac- 
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DP point 
>- (i~1.27) 




FIG. 10: Phase diagram of the system in the (pi,p2) pa- 
rameter space with p su b s = Pc — 0.287338. The dashed line 
(schematic) is a critical line on which quantities in the sys- 
tem are power-law distributed. Above the line there is annular 
growth, and below is sub-critical growth. The bulk exponents 
6, z, v\\ are equal to DP values all along this line, whereas the 
spreading exponents S, rj, r vary continuously along the line 
(representative values of r are indicated). The crosses corre- 
spond to those point along the critical line that were numer- 
ically determined. 



tive site. In this case space-time is divided in two regions: 
the active one, and the outer, inactive one. In Fig. [TT] 
we show a snapshot of a AM growing cluster at a given t. 
We see that sites that make the difference with usual DP 
are mostly located at the boundary of the active region. 
We consider a large box of size £± <C £j_ in space and 
t\i <C £|| in time, sufficiently far away from the boundary 
with the inactive region. Its statistical properties will be 
completely independent of its precise location and are in- 
distinguishable from those of a box with the same size, 
with the fully active initial condition. This shows that 
the bulk exponents 0,vxi v \\ an d thus z are unchanged 
by the compensation process also if we use the single 
seed initial condition. However, the spreading exponents 
5, 77, t are naturally defined only in the seed initial condi- 
tion, and involve averages over all runs. These exponents 
depend continuously on pi and P2 ■ 

A similar scenario happens in 1-dimensional models 
which display critical behavior, despite their breaking of 
the time-reversal symmetry. In these models jjj [101 HH" 
|2"U] each site is active or inactive, as in DP, but is 
equipped with an additional auxiliary field <$> which al- 
lows for a large degeneracy of the absorbing state. We 
discuss DP with auxiliary fields using the example of the 
Threshold Transfer Process (TTP) [TU1 [T7]. In the 1- 
dimensional TTP, a site may be vacant, singly or doubly 
occupied, corresponding respectively to states <7j = 0, 1 
or 2. The auxiliary field <fi denotes the density of singly 
occupied sites. A doubly occupied site corresponds to 
the active state. Contrary to the DP case, there are 
infinitely many absorbing states since any configuration 
with no doubly occupied site is absorbing. Initially, only 
the site at the origin is doubly occupied, while the state 
a, of each other site is 1 with probability 0i n ; t , and oth- 
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FIG. 11: A snapshot of a growing cluster in the AM at com- 
pensation (pi = 0.01, P2 = 0.600). Red (respectively black, 
and blue) points are those sites that have been attempted 
once (respectively twice, and more than twice). 



erwise. At each time step, a site i is selected at random. 
If <7j(i) < 2, then <jj(i + 1) = 1 with probability r and 
with probability 1 — r, irrespective of the precise initial 
value. If <Tj(t) = 2, the site releases one particle to all 
neighbors with a(t) < 2. 

In TTP, r plays the role of control parameter, and in 
d — 1, r c — 0.6894 [T7J. At criticality the bulk exponents 
and the hyperscaling relation behaves as in DP indepen- 
dently of the initial condition. However the spreading 
exponents continuously depend on the initial condition 
(/>i n it. Setting the initial density of singly-occupied sites 
to its stationary value cj> st — r c , one recovers the full 
set of DP critical exponents [10,. As far as we know, a 
theoretical explanation for the continuous change in the 
spreading exponents S, rj is still an open question. 

It is worth mentioning a second class of models with 
similar behavior, which corresponds to DP with special 
absorbing boundary conditions. In particular DP with 
absorbing walls at positions x(t) = ±C • t z shows spread- 
ing exponents that continuously depending on C [51] [55] . 
Analogous results with a moving active wall are presented 
in [22j ■ Moreover, models with soft or hard walls con- 
ditions in one dimension can be studied analytically in 
the case of compact DP, a simplified version of standard 
DP. These models can be mapped onto compact first 
attempt (soft walls) and compact first infection (hard 
walls) . Dickman showed [24] that in this case the critical 
behavior is maintained when pi is reduced, i.e. in this 
case we do not have a stretched exponential behavior. 



In conclusion, memory effects in immunization prob- 
lems, or the presence of auxiliary fields in TTP-like mod- 
els, introduce high degeneracy of the absorbing state and 
thus breaks the time reversal symmetry. In these sys- 
tems, at criticality, the bulk DP exponents arc recovered. 
However, if the initial condition is sufficiently far from its 
stationary value (which is <p st for TTP-like models, and 
the fully twice-attempted lattice for the compensation 
model, or the fully once-infected lattice in the modified 
first Infection Model), the spreading exponents depend 
continuously on the initial condition. Non-stationarity 
seems to play a key role in the observed anomaly of the 
spreading exponents. 

VI. CONCLUSIONS 

We have shown that DP universal behavior is strongly 
affected by changes in the first probabilities to activate 
sites. This modification corresponds to a special case 
of "long memory" , where each site remembers exactly 
how many times it has been activated (or attempted) 
before. Our main result is that, although the change of 
the very first attempt probability takes the model out of 
criticality, by changing the second attempt probability in 
the opposite direction, we can restore critical behavior, 
in a process we call compensation. 

We have concentrated in this paper in the case of two 
spatial dimensions, but qualitatively the same behavior is 
obtained for one dimension. However in one dimension 
the deviations from criticality when p\ is changed are 
much weaker than in two dimensions, then making the 
determination of the compensation condition much more 
difficult numerically 

It is uncertain for us at present if the phenomenology 
of the compensation in the Attempt Model applies also 
to the Infection Model. It seems that the compensation 
effect can occur, but we do not have enough numerical 
evidence to assure that the systematics of critical expo- 
nents is the same as that discussed in this paper for the 
Attempt Model. 

The quantitative change we have obtained of the crit- 
ical exponent r is rather weak. For instance we have 
obtained changes of the value of r of approximately ±0.1 
around the DP value r — 1.268. These variations are typ- 



ically within the error bars of experimentally determined 
exponents in concrete situations. Although changes in 
other exponents like rj and S were found to be larger, in 
some practical situations where avalanches are observed, 
it is the size distribution exponent r the one that can be 
directly measured. So it seems dubious that the effect we 
have discussed can be observed in a concrete realization. 

In this respect we want to mention that we have found 
other realizations of the DP process where the effect is 
quantitatively much more important. For instance, the 
process in which we try to activate neighbors with proba- 
bility p, having also a self-activation probability pg of the 
same site, belong also to the DP universality class. In this 
case we have observed that a lower probability to activate 
neighbors for the first time can also be compensated by 
larger self-activation probabilities during the next steps, 
and in this case the quantitative effect is much more im- 
portant. In particular we have obtained avalanche distri- 
butions with t ~ 1.7, which compares well with those val- 
ues observed for earthquakes. In fact, we have evidence 
that models that have been recently used to justify such 
values of r and other properties of earthquakes [2SJ US] 
can be partially re-interpreted as modified DP models of 
the kind discussed here. We will address this issue in a 
forthcoming presentation. 

Aside from any application to a concrete situation, we 
want to stress the fact that the present model provides 
a link between two classes of models with very different 
behavior. Although we obtained the same results as in 
models with an auxiliary field (criticality with scaling re- 
lations preserved, time- reversal symmetry breaking), our 
microscopical description fits in the framework of mod- 
ified first Infection Models. This may be an interest- 
ing approach to the open problem of initial-condition- 
dependent exponents in absorbing phase transitions. 
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